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I. INTRODUCTION 

Recently, both the singly ionized barium (Ba+) and 
(Ra"*") are proposed as suitable candidates for the parity 
non-conservation (PNC) experiments These candi- 

dates are also suitable for the atomic clock experiments 
[1) 05 Hi- Accurate determination of the electric dipole 
(El) matrix elements are essential in achieving sub-one 
percent PNC amplitudes in the above candidates 0, @ . 
They are also crucial in determining the dipole polariz- 
abilities accurately which are required to estimate shifts 
due to the stray electric fields applied 0] during the 
measurements in the above proposed experiments. There 
is no direct method to measure them experimentally. 
Mainly, these results are estimated from the branching 
ratio (lifetime) measurements from various states. How- 
ever, when branching ratios of two or more different chan- 
nels from a given state are relatively significant, it is not 
possible to estimate the El matrix elements precisely. 
There are a few branching ratio [1, [^, [l3] and lifetime 
prl . [T^ . [ist measurements of the first P-states available 
in Ba"*", although they are not very precise to estimate 
the El matrix elements accurately. On the otherhand, 
such experiments are not available yet in Ra+. 

It is also possible to calculate the above matrix el- 
ements accurately using a highly potential many-body 
method. Ba+ and Ra^ are heavy systems implying both 
the electron correlation and relativistic effects will be 
much stronger in these systems. Variety of relativistic 
many-body methods have already been employed to cal- 
culate the El matrix elements in the considered systems 
Bin, HE [Hi- But the results obtained by aU these 
works differ at significant precision for which they can be 
used to determine sub-one percent accuracy of the PNC 
amplitudes or other experimentally measured quantities 
like light-shift ratios in the considered systems. 



In this work, we have employed the relativistic coupled- 
cluster (RCC) method to account more electron correla- 
tion effects through the non-linear terms for the determi- 
nation of many El matrix elements in both Ba^ and Ra^. 
These results are then used to evaluate the light-shift ra- 
tios in Ba+ at two wavelengths on which experiments 
were conducted [l^, . By analyzing the trends of con- 
tributions to the light-shift ratios and transition proba- 
bilities, we try to give the upper limits to some of the 
important matrix elements by combining with the exper- 
imental result in that system. The matrix elements are 
further used to determine the lifetimes and static dipole 
polarizabilities in Ba+. By following the same procedure, 
we determine in Ra+ the lifetimes of the first excited P- 
states and light-shift ratios at the same wavelengths in 
which experiments were carried out in Ba^. 



II. THEORY AND METHOD OF 
CALCULATIONS 

The shift of energy in an atomic state I7, J, M), where 
J represents the total angular momentum of the state 
with its azimuthal component M and 7 is the additional 
index representing other required quantum numbers, due 
to non-resonant ac light in an average period of light 
oscillations and neglecting the mixing of the magnetic 
sub-levels is given by U 
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where E is the applied vector electric field and is its 
magnitude in the z-component. Here ao, ol\ and are 
the scalar, vector and tensor polarizabilities of the state 
I7, J, M). This expression differs by the term containing 
vector polarizability when the dc electric field is applied. 
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Using the first-order time-dependent perturbation the- 
ory and assuming large detuning from resonance with 
the applied electric field, the scalar, vector and tensor 
dynamic polarizabilities of I7, J, M) given as [2l|, [23j 
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respectively, where D is the El operator, [J] is the de- 
generacy factor equal to 2 J -f 1, 7' and K represent dif- 
ferent states than I7, J, M) with opposite parity, matrix 
elements with double bars represent the reduced matri- 
ces, Ej/i( represent energies of the corresponding states 
in the absence of the external electric field and uj is the 
frequency of the applied electric field. 

The scalar factors (j>i{J, i^)'s (i = 0, 1, 2) multiplied in 
the above equations are defined as 
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The static dipole polarizabilities can be determined by 
substituting w as zero in the above equations. 

In fact, the vector shift is maximal than other two 
shifts for the purely circular polarized light when it is 
aligned with any existing magnetic field [Ij, [IB] ■ In this 
case, it is possible to measure the ratio of light shifts 
instead of measuring the light shifts or light intensities 
directly. Mathematically, it is equivalent to say as the 
following; the ratio of the shift in energies of two different 
states (|7i, Ji,Mi) and I7/, Jj, Mf)) is given by 
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which is known as the light-shift ratio in the literature 

MM- 

The probability coefficient (in s^^) due to the allowed 



transition is given by 
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where A (A) and Sflii= \{f\\D\\i}\^) (au) are the wave- 
lengths and strengths of the corresponding transitions, 
respectively. 
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FIG. 1; Schematic low-lying energy level diagrams and decay 
of the P-states in Ba"'' and Ra^. n = 6 and n = 7 in Ba^ and 
Ra"'", respectively. 

As seen in Fig. [1] the first excited P-states in both Ba"*" 
and Ra+ can decay either to the ground S-states directly 
or via the excited D-states through the allowed transi- 
tions. The lifetimes of these P-states can be evaluated 
from the inverse of the corresponding net transition prob- 
abilities due to all possible transition channels. There- 
fore, the net transition probabilities for the P-states are 
given by 



- ^nPi/2-»nSi/2 + ^nPi/2^("-l)D3/2(2-10) 
^"P3/2 = ^nP3/2— nSi/2 + ^nP3/2^(n-l)D3/2 

+^"P3/2-(n-l)D5/2 (2-11) 

and hence, the lifetimes (in s) of these states are given 
by 
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where n — 6 and n = 7 in Ba+ and Ra"*", respectively. 

Since the transition wavelengths can be determined 
from the excitation energies from the corresponding tran- 
sitions, it is clear from the above equations that both the 
polarizabilities and transition probabilities depend upon 
the excitation energies (or wavelengths) and El matrix 
elements. Since our motivation in the present work is to 
verify the accuracies of the El matrix elements, we con- 
sider the experimental (observational) excitation energies 
(wavelengths) for our purpose and analyze only the role 
of the El matrix elements from our calculations. 

We employ the RCC theory for the Dirac-Coulomb 
Hamiltonian in the single and double excitations ap- 
proximation along with the important triple excitations 
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(CCSD(T) method) to calculate the atomic wave func- 
tions in both Ba^ and Ra"*". In this procedure, we express 
the one valence (v) configuration state function as [2g.[27j 

= e^{l + 54|$,), (2.14) 

where the reference state is constructed by append- 
ing the corresponding valence orbital to the closed-shell 
Dirac-Fock (DF) wave function |$o); which in second 
quantization formalism is expressed as \^v) — o-t\^o)- In 
the above equation, T and Sy are the core and valence- 
core correlation operators. In the linear approximation. 



the above equation yields the form 

= {1 + T + 5,„}|$,). (2.15) 

As shown in the earlier works, contributions (electron 
correlation effects) from higher triple and quadrupole 
excitations arise through the non-linear terms in the 
CCSD(T) method and their contributions to the various 
properties in the heavy systems are significantly large 

The reduced transition matrix elements 
between two states jvpy) and are evaluated by 
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where D — e^^ De^ and N — e^ 'e^ are the non- 
truncative series in the ordinary RCC theory. However, 
in the CCSD(T) method, only up to the effective three- 
body terms, when they are expanded using the general 
Wick's theorem, of these series will contribute to the 
above expression. These terms are calculated step by 
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step and sandwiched between the necessary Sy operators 
to evaluate the matrix elements. In order to estimate the 
errors associate with the calculated El matrix elements, 
we check the convergence of these results using both the 
length and velocity gauge expressions which at the single 
orbital level are given by 
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in the velocity form, respectively. In these expressions, 
C*^^-* is the Racah vector with rank one, Pi/f{r) and 
Qi/f{r) are the large and small components of the Dirac 
single particle orbitals, Ki/f are the relativistic quantum 
numbers, Ci/f are the single particle energies and a is the 
fine structure constant. 



III. RESULTS AND DISCUSSIONS 
A. Ba+ 

We present our calculated El matrix elements of Ba"*" 
in TablelH The reported error bars in these results except 
from the matrix elements related to the D- and F- states 
are estimated from the differences between the length and 
velocity gauge results. It was rather difficult to converge 



the velocity gauge results for the D- and F- states, there- 
fore we consider the difference between the converged 
length gauge results using the CCSD(T) method and just 
by single and double excitations approximation (CCSD 
method), which are assumed as the upper limit to the 
contributions due to the neglected excitations, as the er- 
ror bars. However, these error bars are just used to es- 
timate errors in the light-shift ratio calculations in the 
considered system and we analyze the accuracy of these 
results later by fitting with the experimentally measured 
values of the light-shift ratios at two different wavelengths 
as discussed below. 

We also compare our results with the earlier works in 
Table IT] As seen, most of the important low-lying El 
matrix elements from various works differ at the second 
decimal places which means the difference is mainly due 
to the different many-body methods employed and they 
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TABLE I: Absolute magnitudes of the reduced dipole matrix 
elements in Ba+. Estimated error bars from this work given 
inside the parenthesis. 



Transition Present Ref. [14.1 Ref- [15] Ref. [16] Ref. [VT} 



QPl/2 


-» 


6Sl/2 


3.36(1) 


3.310 


3.3266 


3.3357 


3.300 


7Pl/2 




6Sl/2 


0.10(1) 


0.099 


0.1193 


0.0621 




8pi/2 




6Si/2 


0.11(5) 


0.115 


0.4696 






6P3/2 




6Sl/2 


4.73(3) 


4.674 


4.6982 


4.7065 


4.658 


7P3/2 


— > 


6Sl/2 


0.17(5) 


0.035 


0.3610 


0.0868 




8p3/2 


— > 


6Sl/2 


0.11(5) 


0.073 


0.5710 






6pi/2 




7Sl/2 


2.44(4) 


2.493 


2.3220 






6pi/2 




8Si/2 


0.66(5) 


0.705 


0.7283 






6P3/2 


-> 


7Si/2 


3.80(2) 


3.882 


3.6482 






6P3/2 




8Sl/2 


0.97(5) 


1.025 


1.0518 






6pi/2 




5^3/2 


3.11(3) 


3.055 


2.9449 




3.009 


7pi/2 


-> 


5^3/2 


0.28(2) 


0.261 


0.3050 






8pi/2 




5^3/2 


0.13(2) 


0.119 


0.1121 






6P3/2 


— > 


5^3/2 


1.34(2) 


1.334 


1.2836 




1.312 


7P3/2 




5^3/2 


0.16(1) 


1.472 


0.1645 






8p3/2 




5^3/2 


0.07(2) 


0.070 


0.0650 






4/5/2 




5^3/2 


3.75(11) 










5/5/2 




5^3/2 


1.59(8) 










6/5/2 




5^3/2 


0.17(2) 










6P3/2 




5^5/2 


4.02(7) 


4.118 


3.9876 




4.057 


7p3/2 




5^5/2 


0.46(1) 


0.432 


0.4788 






8P3/2 




005/2 


0.21(2) 


0.206 


0.1926 






4/5/2 




5^5/2 


1.08(4) 










5/5/2 




5^5/2 


0.45(7) 










6/5/2 




5d5/2 


0.15(2) 










4/7/2 




5^5/2 


4.84(5) 










5/7/2 




5^5/2 


2.47(6) 










6/7/2 




5d5/2 


1.04(7) 










6pi/2 




6^3/2 


4.89(10) 










6pi/2 




7(^3/2 


1.50(8) 










6P3/2 




6^3/2 


2.33(7) 










6P3/2 




7^3/2 


0.67(4) 










6P3/2 




6(^5/2 


6.91(21) 










6p3/2 




7(^5/2 


2.01(5) 











may not be due to the numerical methods used in the 
calculations. That is the reason why we have presented 
our results up to the second decimal places in the above 
table. As a matter of fact, correct results up to the sec- 
ond decimal places are crucial in obtaining the light-shift 
ratios correctly as will be shown below. First, we would 
like to briefly discuss here the important differences in 
the various works presented in Table [J Geetha et al. [isj 
have employed the same approach like ours, but they had 
used part numerical orbitals using GRASP and part ana- 
lytical orbitals using Gaussian type orbitals (GTOs). In 
the present work, we have considered pure analytical or- 
bitals using GTOs. Another difference between these two 
works is that D of Eq. (|2.16p is truncated only at the 
effective one-body terms by Geetha et al. [3^, where 
other higher order terms have been considered in the 



present calculation. The main difference in the methods 
applied by Dzuba et al. [l3] and ours have been discussed 
by Geetha et al. in their work [l^. In summary, they 
have employed the Green's function technique which is 
also an all order perturbative method like ours. In their 
work, the orbitals are also obtained analytically like the 
present work. However, the level of approximation in 
both the works could be different which cannot be obvi- 
ously distinguished. Iskrenova-Tchoukova and Safronova 
[T^ have used linearized CCSD method with important 
triple excitations and orbitals have been constructed us- 
ing the B-spline basis. As has been pointed out earlier, 
the non-linear terms can incorporat e hig her order correla- 
tion effects which are necessary [H, [29f to account in the 
considered heavy system. Again, the procedure to con- 
sider partial triple excitations in their method and in the 
present work is different. We consider effects of partial 
triple excitations by estimating each time their contribu- 
tions to the energy of the corresponding valence state and 
then solve the CCSD method amplitudes self-consistently 
which involves the above energy. However, Iskrenova- 
Tchoukova and Safronova include them directly in the 
amplitude determining equations. Guet and Johnson 
have just applied the relativistic many-body method at 
the second order approximation (MBPT(2) method) us- 
ing the B-spline basis to calculate their results in contrast 
to our method which is all order in nature [l^, . 

We have used our calculated reduced El matrix el- 
ements and experimental energies [3l| in Eq. (|2.3p to 
obtain the dynamic vector dipole polarizabilities (ai) of 
the 6S1/2 and 5D3/2 states in Ba^ at two different wave- 
lengths (of external electric fields) and they are presented 
in Table M Only 6P1/2, 7Pi/2, 8P1/2, 6P3/2, 7P3/2, 
8P3/21 4F5/2, 5F5/2 and 6F5/2 intermediate states (K) 
have been taken into account as they seem to be the 
most important contributing states as seen from this ta- 
ble. It can be noticed that the magnitude of ais of the 
6S1/2 state are larger than the 5D3/2 state at both the 
wavelengths, but their differences are smaller when the 
frequency of the external field is small. Again, the sign 
of these quantities are different at the smaller wavelength 
(514.53 nm), whereas they are opposite at the large wave- 
length (1111.68 nm). For both the states, the largest 
contribution comes from the 6P1/2 state. The next high- 
est contribution arises from the 6P3/2 state. The third 
largest contributions to ai of the 5D3/2 state comes from 
the 4F5/2 state. Surprisingly, the rest of the contribu- 
tions are dominated by the higher F-states compared to 
the P-states. We have given contributions from pure core 
correlations, core-valence correlations and higher states 
those are neglected as " Others" . The core correlation ef- 
fects are around 0.319 au and 0.147 au at the wavelengths 
514.53 nm and 1111.68 nm, respectively. These contri- 
butions are evaluated using the MBPT(2) method. The 
possible errors due to the omission of the higher states 
may come from the neglected F-states which seems to 
be important in the determination of ai of the 5D3/2 
state. Therefore, an approach similar to the one applied 
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to determine the dynamic scalar and tensor polarizabil- 
ities [m will be more appropriate to use for the high 
precision light-shift ratio studies in the present system. 

Now considering the above results of ai for both 
the 6S1/2 and 5D3/2 states and using Eq. (|2.8p . we 
obtain the hght-shift ratios in Ba+ as i?A=5i4.53nm ^ 
-11.36(20) and i?>=iiii-68nm ^ o.418(3), where core cor- 
relations contribute around 0.6% and 1% at the wave- 
lengths 514.53 nm and 1111.68 nm, respectively. The 
corresponding experimental results are Jj-*'=5i4.53nm _ 
-11.494(13) and i?A=iiii.68nm ^ o.4176(8) at the wave- 
lengths 514.53 nm and 1111.68 nm, respectively p^[20l|. 
Our results match well with the corresponding experi- 
mental results. It shows that accounting core correla- 
tion effects are essential for the high precision results in 
the light-shift studies, although their magnitudes seem 
to be very small in comparison to the ai of the 6S1/2 
state. Especially, their contributions are relatively larger 
at the higher wavelength (smaller frequency) . The errors 
quoted in our light-shift ratio calculations are obtained 
from the error bars given for the El matrix elements in 
Table D 

We are now in a position to analyze the accuracy of 
the El matrix elements by studying their contributions 
to the light-shift ratio calculations. Since contributions 
from the 6P1/2 state seems larger, we first analyze ma- 
trix elements involving this state. From Table U we can 
find that different theoretical predictions of the magni- 
tude (sign is irrelevant for our observable) (6pi/2||-D||6s) 
range from 3.30 au to 3.37 au (along with the error bar). 
Now keeping all other matrix elements unchanged from 
our calculations, we observe that light-shift ratios vary 

from i?^=514-53nm ^ _10.77and i?A=1111.68nm ^ 0.337 to 
^A=514.53nm ^ -HAh and i?A=1111.68nm ^ 0.432. With 

(62i/2||-D||6s) = 3.33, which is reported by Geetha et al. 
[isl l and Iskrenova-Tchoukova and Safronova [16^, we ob- 
serve i?^=514.53«m ^ _11.06 and i?A=llll-68nm ^ 0.337. 

Keeping the value of the above matrix element con- 
stant, we consider now the range of (6p3/2||^||6s) as 
4.66 au to 4.76 au and observe from i?^=5i4.53nm ^ 

-11.21 and i?A=1111.68«m ^ 0.432 to i?^=514.53nm ^ 

-10.99 and i?^=iiii-68nm ^ Q 353^ improve these 
results, we fix now (6p3/2||^||6s) = 4.70 and vary 
(6pi/2||-D||6s) result again. We find that when we con- 
sider (6pi/2||i:'||6s) = 3.35 and (6pi/2||i:'||6s) = 3.37 for 
(6P3/2||-D||6s) = 4.70, we get i?^=5i4.53nm ^ -11.32 and 

^A=1111.68nm ^ q 433 ^nd i?A=514.53nm ^ _11.51 and 

^A=iiii.68nm ^ 0.455, respectively. Since the former 
result, is little close with the experimental results for 
both the wavelengths, we now keep (6pi/2||-D||6s) = 3.35 
and increase (6p3/2||i?||6s) result and get i?A=5i4.53nm ^ 

-11.27 and i?A=1111.68nm ^ 0.412 and i?^=514.53nm ^ 

-11.22 and i?A=iiii.68«m ^ o.396 for 4.72 au and 
4.74 au, respectively. Therefore, we conclude now that 
(6p3/2||-D||6s) will be around 4.73 au. Now keeping 
(6pi/2||-D||6s) = 3.35 and {Qp^/2\\D\\&s) = 4.73 when 
we vary (6pi/2||-D||5c?3/2) from 3.00 au to 3.15 au, we 



get from i?A=514.53nm ^ -11.93 and i?A=1111.68«m ^ 
0.435 to i?^=514.53«m ^ -11.02 and i?>=llll-68nm ^ 

0.394. Results close to the experimental measure- 
ments are able to produce only when (6pi/2||^|| 5^3/2) 
is around 3.08 au. Now we keep (6^1/2! I-DII6S) = 3.35, 
(6]53/2||-Ci||6s) = 4.73 and (6pi/2||i:'||5(i3/2) = 3.08 then 
vary (6^3/2 5^3/2) matrix element from 1.28 au to 
1.36 au. We get light-shift ratios close to experimental 
results when (6^3/2 ||-D|| 5^3/2) is around 1.32 au. How- 
ever, light-shift results vary slowly in the given range 
of this matrix element while the previous matrix ele- 
ments considered earlier play the crucial roles in de- 
ciding the final results. We now vary (4/5/2 ll^*! 15^3/2) 
matrix element from 3.50 au to 4.00 au by consid- 
ering (6^1/2! |Z?||6s) = 3.35, (6p3/2||^||6s) - 4.73, 



(6p 



1/2 



\D\\5d: 



3/2/ 

we get from i?A=514.53mn ^ -11.57 and i?A=1111.68nm ^ 
0.409 to i?^=514.53«m ^ -11.36 and i?A=1111.68nm ^ 

0.418. By fixing (4/5/2 1 1-0115^3/2) = 4.00, when we 
reshuffle other matrix elements, we get light-shift ra- 
tios close to the experimental values for (6pi/9||Z?||6s) = 
3.35, (6p3/2||Z?||6s) = 4.72, {6pi/2\\D\\5dy2) = 3.08 
and (6p3/2 1 l-Dj 15^3/2) = 1.34. When we keep these ma- 
trix elements constant and vary the matrix element of 
(4/5/211^115^3/2), we get the best light-shift ratios for 
3.65 au and they correspond to i?^=5i4-53nm ^ -11.49 
and ij^=iiii-68nm ^ 0.419. For any other combinations, 
the light-shift ratios vary by large amount from the ex- 
perimental results; at least for one of the wavelengths. 

Again, there are also measurements of the transition 
probabilities in Ba"*" available from three different exper- 
iments [1, We evaluate them using the experi- 
mental wavelengths (determined from the corresponding 
experimental excitation energies jsij) and El matrix el- 
ements presented in Table [D These results are given in 
Table IIIII We have also presented both the experimen- 
tal and other theoretical results in the same table. It 
can be clearly noticed that our results in all the cases 
match well with the experimental results with smaller 
error bars where as some of the earlier works differ sig- 
nificantly. Differences between all the theoretical works 
are discussed by Geetha et al. [Tsj . 

Now if we consider (6pi/2||i5||6s) as 3.35 au as 
we analyzed from the light-shift ratio studies, we get 
746p^/2^6Si/2 94.566 x 10^s~^ which is at the marginal 
lower side error limit of Reader et al experimental re- 
sult [9]. Our ^gPj^2^5D3/2 result is also at the marginal 
upper limit of Kastberg et al experimental result [1] . If 
we consider {6pi/2||-D|| 5^3/2) = 3.08 from the analysis 
of light-shift ratio studies then we get ^6Pi/2^5D3/2 = 
35.015 X lO^s""'^ which is within the experimental error 
bar. 

By gathering all the informations from the analysis of 
roles of various matrix elements in the light-shift ratios 
and transition probabilities studies, we arrive in the con- 
clusion at this point that the magnitude of the following 



3.08 and (6^3/2 ||L'||5d 



^3/2/ 



1.32, 
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matrix elements will be 



(6pi/2||i^||6si/2) 
(6p3/2l|i5||6si/2) 
(6pi/2||D||5d3/2) 

(6p3/2 1 1-01154/2} 

(4/5/2||D||5d3/2) 



3.36(2) 
4.73(3) 
3.08(3) 
1.34(2) 
3.73(20). 



The above conclusions differ from the findings of the pre- 
vious work by Sherman et al. [2d| who find 

(6pi/2||i^||5d3/2) = 3.14(3) 
(4/5/2||7^||5d3/2) = 4.36(36), 

where they had combined with Geetha et al [l^ results 
to estimate these values. Further more, light-shift ratio 
measurements with a large number of wavelengths are re- 
quired in order to estimate the above El matrix elements 
more accurately. 

There are a number of measurements on the lifetimes 
of the 6p ^Pi /9 and 6p ^P3/2 states are available in the 
literature [Til. Tl^. ITsj . but the associate error bars of 
these measurements are very large. We evaluate the net 
transition probabilities of these states using the above 
matrix elements and experimental wavelengths [3l| . We 
give these results in Table IIVI and compare with the cor- 
responding experimental results. Results obtained from 
other works are also given in the same table. As it can be 
seen from this table that the lifetime of the 6P1/2 state is 
more than 2% difference from the earlier findings which 
needs to be verified by new experiments. Although, our 
calculation of the lifetime of the 6P3/2 state is within 
the error bar of the experimental results and agree with 
other findings, but the Agpg^j ^505/2 contributes around 
28% in its evaluation. Since the accuracy of the El ma- 
trix element of (6^3/2 l|i?||5d5/2) is not verified well as 
it has been done for the transition matrix elements in- 
volved in the determination of lifetime of the 6P1/2 state, 
the lifetime of the 6P3/2 state may be little different than 
what we find, but it will be within the error bar that we 
have given (which is large). The error bars given in our 
lifetime calculations come only from the error bars of the 
El matrix elements as given above. 

Since the knowledge of the static dipole polarizabili- 
ties of different states in Ba"*" is important to estimate 
shifts in the energy levels when they are subjected to 
high precision experiments like atomic clock, PNC and 
so on, we have calculated these quantities for the first 
five low-lying states in the given system. Likewise vec- 
tor polarizabilities, the scalar and tensor static dynamic 
polarizabilities are evaluated using our El matrix ele- 
ments presented in Table |T] and their experimental exci- 
tation energies [Ui]. We have considered the 6S1/2, 7Si/2, 
8S1/2, 5D3/2, 6D3/2, 7D3/2, 5D5/2, 6D5/2, 7D5/2, 4F7/2, 
5F7/2 and 6F7/2 states in these calculations along with 
the intermediate states considered for the above vector 
polarizabilities calculations. Again, we have calculated 
the core correlation effects using the relativistic CCSD 



method employed for the closed-shell system in our ear- 
lier work |33i]. The core- valence and higher state contri- 
butions other than the matrix elements reported in Table 
Uare evaluated using the MBPT(2) method. Our static 
dipole polarizability results of various states in Ba+ are 
presented in Table |Vl The core-correlation effects for 
the scalar and tensor polarizabilities are 9.582 an and 
—0.372 au, respectively. There are experimental results 
available only for the ground state [H, [s^. Again, a 
couple of theoretical calculations available for the ground 
state [H, [H, [13, llll . Except Iskrenova-Tchoukova and 
Safronova work [16j |. others have employed molecular 
codes to determine them. Also, Iskrenova-Tchoukova and 
Safronova have used the El matrix elements from the lin- 
earized CCSD(T) method to evaluate the ground state 
polarizabilities. In their work, they have considered only 
four important matrix elements to evaluate them and we 
have considered two more higher excited states matrix el- 
ements in this calculation. Again, their core-correlation 
effects are considered using the lower order many-body 
methods and it is larger than our finding. There are no 
other results for the excited states available to compare 
with our results. The error bars given in our calculations 
are from the errors given for the El matrix elements in 
Table H] The agreement between our ground state polar- 
izability result with its experimental measurement fur- 
ther supports the accuracy of the El matrix elements 
used in this calculation. 



B. Ra+ 

We have also followed the same procedure as above to 
determine the light-shift ratios and lifetimes of the first 
excited P-states in Ra+. We present the various El ma- 
trix elements used in these calculations in Table IVTl The 
error bars of different matrix elements are estimated us- 
ing the same procedure followed as in Ba"*". There are also 
other calculations available in this system Dzuba 
et al. [IJ have used the same method as in Ba"*" to cal- 
culate these matrix elements. The procedure followed by 
Safronova et al. |l8i] is based on the linearized CCSD(T) 
method like Iskrenova-Tchoukova and Safronova work in 
Ba"*" jJJl]. Most of our results match with these calcula- 
tions, but they still differ at the second decimal places as 
in Ba+. 

Since Ra"*" has been proposed for the atomic clock Q 
and PNC experiments 0, knowledge of the accuracies 
of the above El matrix elements are essential. We pro- 
pose also similar light-shift ratio measurements at differ- 
ent wavelengths for the same purpose. As a preliminary 
study, we consider the same off-resonant wavelengths as 
have been considered in the Ba+ experiments [1^ i2^ 
and have carried out light-shift ratio calculations. In Ta- 
bic IVII[ we present the dynamic vector polarizabilities 
for the 7Si/2 and 6D3/2 states at A = 514.53nm and 



A = 1111.68nm. We use our El matrix elements pre- 
sented in Table IVII with the experimental energies [3l| 
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to evaluate them. In this calculation, we have consid- 
ered the 7Pi/2, 8P1/2, 9Pi/2, 7P3/2, 8P3/2, 9P3/2, 5F5/2, 
6F5/2 and 7F5/2 intermediate states. 

As seen in Table IVIIl the trend of the contributions 
from various intermediate states is same as in Ba"*". The 
core-correlations are little larger compared to Ba+ and 
they are around 0.346 au and 0.159 au at A = 514.53nm 
and A = 1111.68nm, respectively. We obtain around 

^A=514.53mn ^ -6.42(7) and i?A=llll-68nm ^ 0.017(1). 

The given error bars again come from our estimated ac- 
curacies of the El matrix elements. Since there are no 
experimental data available in Ra"*" to check the accuracy 
of the El matrix elements or light-shift ratios, our results 
will be useful for the future light-shift experiments in this 
system. 

We have already reported polarizabilities of the 7Si/2, 
6D3/2 and 6D5/2 states in our earlier work 5]. How- 
ever, the knowledge of the transition probabilities and 
lifetimes of the first P-states in the considered states are 
not known yet. We use our El matrix elements reported 
in Table IVII and experimental energies [3l| to evaluate 
them. In Table IVlIIl we present the transition probabili- 
ties of various channels from the 7Pi/2 and 7P3/2 states. 
As seen from this table, the relative transition probability 
of the 7P3/2 7Si/2 transition is much larger than the 
7Pi/2 ^ 7Si/2 transition in the present system compared 
to Ba+. Since the wavelength of the 7Pi/2 6D3/2 tran- 
sition is one order larger here, the corresponding tran- 
sition probability is comparatively smaller. Hence, the 
lifetime of the 7Pi/2 state becomes little larger than the 
6P1/2 state in Ba+. Again, the lifetime of the 7P3/2 
state becomes smaller than the 6P3/2 state in Ba+ due 



to the fact that the 7P3/2 7Si/2 transition probability 
is much larger. We present the net transition probabili- 
ties and lifetimes of the 7Pi/2 and 7P3/2 states in Table 
HXl In Ra+, A7P3^2^6D5/2 contributes only by 12% to the 
final lifetime determination of the 7P3/2 state in contrast 
to 28% in Ba+. 



IV. CONCLUSION 

We have employed the relativistic coupled-cluster the- 
ory to calculate the electric dipole matrix elements in 
the singly ionized barium and radium. These elements 
are used to determine the light-shift ratios and transition 
probabilities in the singly ionized barium and comparing 
with the corresponding experimental data, we have esti- 
mated the accuracy of various low-lying matrix elements 
in this system. Further, we have evaluated the scalar and 
tensor static dipole polarizabilities and lifetimes in the 
same system. We have also calculated the light-shift ra- 
tios, transition probabilities and lifetimes using our elec- 
tric dipole matrix elements in the singly ionized radium. 
These data will be helpful in the future experiments in 
the considered system. 
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TABLE II: Dynamic vector polarizabilities (in au) at two different wavelengths in Ba"*". 

Intermediate 6si/2 5d3/2 

States (K) X = 514.53nm A = 1111. 68nm A = 514.53nm A = 1111. 68nm 



6pi/2 


-978.638 


-45.080 


48.786 


-20.419 


7Pl/2 


-0.014 


-0.006 


-0.035 


-0.014 


8pi/2 


-0.010 


-0.004 


-0.004 


-0.002 


6p3/2 


305.338 


36.717 


6.112 


-1.121 


7P3/2 


0.019 


0.008 


-0.004 


-0.002 


8p3/2 


0.005 


0.002 


-5 X 10"* 


-4 X 10" 


4/5/2 






3.987 


1.541 


5/5/2 






0.153 


0.162 


6/5/2 






0.004 


0.002 


Others 


0.301 


0.137 


0.243 


0.177 



Total -672.994(1.944) -8.223(201) 59.298(1.375) -19.675(323) 



TABLE III: Transition strengths (au), wavelengths (A) and 
probabilities (xlO^s^^) from P-states in Ba+. 

Transition(/ -> i) Sfjj X]^^ Af^^ Others Expt. 



6P1/2 ^ 6S1/2 11.290 4935.5 95.131 93.68° 95(9)^= 

91.78'' 95.5(10)'* 
95(7)" 

22.373 4555.3 119.889 119.37'^ 106(9)'= 
116.25'' 117(4)'* 
118(8)'= 

A6Pi/2--,5D3/2 9.672 6498.7 35.701 32.609'' 33.8(19)'= 

33.42'' 33(8)'' 
33(4)"= 

^6P3/2-,5D3/2 1.796 5855.3 4.531 4.255" 4.69(29)'= 

3.342'' 4.8(5)'' 
4.48(6)'= 

A6P3/2^5D5/2 16.160 6143.4 35.305 34.93'^ 37.7(24)'= 

35.95'' 37(4)'' 
37(4)'= 

^As are determined from the exp erimental excitation 
energies \3B . 
References: "[IE], ''[131, "g], ''[1 and 
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TABLE IV: Net transition probabilities (xlO'^s^^) and life- 
times (ns) in Ba^. 





A ■ 




Expt 










6P1/2 


130.146 


7.68(7) 7.89" 


7.92(8)" 














7.92'' 


7.90(10)-^ 














7.83^= 
















7.99'* 












6P3/2 


159.72 


6.26(11) 6.30" 


6.32(10)-'^ 














6.31'' 


6.312(16)3 














6.27" 
















6.39'* 












References: 






1 and 














TABLE V: Scalar and tensor static {u 


— 0) polarizabilities in Ba+. 






Intermediate 


6S1/2 


6P1/2 


6P3/2 


5D 


3/2 


5D5 


/2 


States (K) 


Scalar 


Scalar 


Scalar Tensor 


Scalar 


Tensor 


Scalar 


Tensor 


6S1/2 




-40.763 


-37.280 37.280 










7Si/2 




19.714 


25.889 -25.889 










8S1/2 




0.844 


0.954 -0.954 










6P1/2 


40.763 






22.992 


-22.992 






7Pl/2 


0.015 






0.064 


-0.064 






8P1/2 


0.014 






0.011 


-0.011 






6P3/2 


74.560 






3.846 


3.077 


24.210 


-24.210 


7P3/2 


0.043 






0.021 


0.017 


0.116 


-0.116 


8P3/2 


0.014 






0.003 


0.003 


0.019 


-0.019 


5D3/2 




-45.984 


-3.846 -3.077 










6D3/2 




68.101 


8.275 6.620 










7D3/2 




4.163 


0.434 0.347 










5D5/2 






-36.315 7.263 










6D5/2 






72.167 -14.433 










7D5/2 






3.895 -0.779 










4F5/2 








11.857 


-2.371 


0.668 


0.763 


5F5/2 








1.643 


-0.329 


0.089 


0.102 


6F5/2 








0.017 


-0.003 


0.009 


0.010 


4F7/2 












13.344 


-4.766 


5F7/2 












2.655 


-0.948 


6F7/2 












0.430 


-0.154 


Others 


9.137 


11.304 


11.716 -0.508 


8.981 


-.019 


9.000 


-.084 


Total 


124.546(1.256) 17.379(2.920) 45.890(3.831) 5.870(128) 49.435(1.458) 


-22.694(538) 50.540(1.410) - 


-29.422(931) 


Expt. [34] 


123.88(5) 
















125.5(10) 














Theory [16] 


124.15 














[36] 


123.07 














[37] 


126.2 














[38] 


124.7 
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TABLE VI: Absolute magnitudes of the reduced dipole matrix 
elements in Ra+. Estimated error bars from this work given 
inside the parenthesis. 



Transition Present Ref. [14] Ref. [18] 



7Pl/2 


-> 7Si/2 


3.28(2) 


3.224 


3.2545 


8Pl/2 


-> 7Si/2 


0.08(4) 


0.088 




9Pl/2 


^ 7Si/2 


0.09(3) 


0.116 




7p3/2 


7Si/2 


4.54(2) 


4.477 


4.5106 


8p3/2 


-> 7Si/2 


0.49(2) 


0.339 




9p3/2 


7Si/2 


0.30(2) 


0.095 




7pi/2 


6^3/2 


3.62(5) 


3.550 


3.5659 


8pi/2 


6^3/2 


0.06(2) 


0.013 




9pi/2 


^ 6^3/2 


0.02(1) 


0.013 




7p3/2 


6^3/2 


1.54(2) 


1.504 


1.5117 


8p3/2 


^ 6^3/2 


0.15(2) 


0.127 




9p3/2 


6^3/2 


0.07(2) 


0.057 




5/5/2 


^ 6d3/2 


4.67(2) 




4.4491 


6/5/2 


^ 6^3/2 


0.86(4) 






7/5/2 


Gds/2 


0.48(11) 






7P3/2 


6^5/2 


4.83(8) 


4.816 


4.8232 



TABLE VII: Dynamic vector polarizabilities (in au) at two different wavelengths in Ra^. 

Intermediate 7si/2 6^3/2 

States (K) A = 514.53nm A = 1111. 68nm A = 514.53nm A = 1111. 68nm 



7pi/2 


-391.483 


-37.764 


31.921 


-869.894 


8pi/2 


-0.008 


-0.003 


-0.002 


-8 X 10"" 


9pi/2 


-0.007 


-0.003 


-1 X 10"* 


-5 X lO"'^ 


7p3/2 


94.788 


22.384 


3.784 


-2.632 


8P3/2 


0.144 


0.059 


-0.005 


-0.002 


9p3/2 


0.035 


0.015 


-6 X 10"* 


-3 X 10"" 


5/5/2 






9.452 


3.361 


6/5/2 






0.149 


0.061 


7/5/2 






0.035 


0.015 


Others 


0.323 


0.145 


0.795 


0.345 



Total -296.207(4.909) -15.167(255) 46.128(679) -868.749(24.027) 



TABLE VIII: Transition strengths (au) , wavelengths (A) and 
probabilities (xlO^s"^) from P-states in Ra^. 

Transition(/ ^ i) gfi^ A/_>i Afj^ 



^7P3/2^6D» ,0 



11.086 


4683.6 


106.083 


20.612 


3815.5 


187.960 


13.104 


10791.2 


10.564 


2.372 


7080.0 


3.385 


23.329 


8022.0 


22.890 



12 



TABLE IX: Net transition probabilities (xlO^s"^) and life- 
times (ns) in Ra"*". 



State(= i) 


Ai 


n 


7Pl/2 


116.647 


8.57(12) 


7P3/2 


214.235 


4.67(5) 



